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Abstract. The radial distribution of globular clusters in 
our Galaxy, M31 and M87 is studied and compared with 
that of halo stars. The globular cluster distributions seem 
significantly flatter than those of the stars bulge. Assum¬ 
ing this is a consequence of an evolution of the globular 
cluster distribution in these galaxies, a comparison with 
the (unevolved) stellar distribution allows us to obtain 
estimates of the number and total mass of clusters lost, 
which are possibly gone to feed the massive central objects 
present in these galaxies. 

It results that the cluster systems studied should have 
been initially about one third and one forth richer than 
now in our Galaxy and in M31, respectively, and twice as 
abundant in M87. The estimated mass of globular clusters 
lost is compatible with the galactic nucleus masses. 
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1. Introduction 

The globular cluster systems (GGS) in the two Virgo gi¬ 
ant galaxies M87 and M49 are clearly less concentrated 
than the halo star distribution (Harris, 1986). Probably 
this feature is not common to all galaxies, even though 
in many cases the available data are probably not good 
enough to compare reliably the cluster and halo distribu¬ 
tions. Anyway, a safe statement is that no case has been 
found where the GCS is more centrally concentrated than 
the halo (Harris, 1991) This is confirmed by the recent 
HST WFPC2 observation of 14 elliptical galaxies (Forbes 
et al 1996). 

A possible explanation of this difference in the distri¬ 
butions has been suggested by Harris & Racine (1979) and 
by Racine (1991) as a difference in the formation ages of 
halo stars and globular clusters. Following these authors. 
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globular clusters are formed earlier, when the density dis¬ 
tribution was less peaked. This possibility cannot be ruled 
out, however it is not supported by any evidence of a sig¬ 
nificant older age for globular clusters respect to the halo: 
this age difference should be large enough to have allowed 
the mother galaxy to contract enough to form halo stars 
in a distribution as more concentrated than globulars as 
observed. Note that in disk galaxies the epoch of cluster 
formation could be early enough to force chemical enrich¬ 
ment but not to take on a distinct spatial structure (Har¬ 
ris 1986). Moreover, this picture does not explain why the 
tails of the two density distributions are about the same. 
Probably, a simpler explanation, working in the majority 
of cases, is the coeval birth of globular clusters and halo 
stars with a further evolution of the GCS radial distribu¬ 
tion, while the collisionless halo stands almost unchanged. 
The causes of evolution are dynamical friction and tidal 
interaction with a compat nucleus; these phenomena can 
act to deplete the GCS in the denser inner galactic re¬ 
gions so to modify the initial radial distribution just in the 
central region leaving unchanged the outer profile, which 
remains similar to that of the halo component. If this is 
true, the halo radial profile clearly represents the shape of 
the initial cluster distribution. 

2. Globular cluster system evolution 

There are various indications that the GCS in a galaxy 
does not behave like a dissipationless system, as the halo 
component is. Actually even if the galaxy is not a spiral 
where disk shocking is an important cause of evolution, 
the tidal shock due to the passage near to the galaxy cen¬ 
tre, and dynamical friction (which acts so to bring mas¬ 
sive clusters closer and closer to the centre) are relevant 
causes of GCS evolution. This leads to a more or less im¬ 
portant change of the GCS spatial distribution and mass 
function. The relevance of the mentioned phenomena de¬ 
pends on the galaxy characteristics: triaxiality enhances 
the efhciences of both of them (Ostriker, Binney & Saha, 
1989; Pesce, Capuzzo-Dolcetta & Vietri, 1992; Capuzzo- 
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Dolcetta 1993 (hereafter CD); Capuzzo-Dolcetta 1996). 
In particular, Pesce et al. (1992) showed that clusters on 
box orbits in a triaxial potential lose their orbital energy 
at a rate one order of magnitude larger than on loop or¬ 
bits of comparable size and energy (even if they are quite 
elongated). It is very likely that a large fraction, if not all, 
of globular clusters are actually moving on box orbits and 
certainly not on quasi-circular loops, due to their early 
formation during the almost radial proto-galaxy collapse 
(see Binney 1988). So previous evaluations of dynamical 
friction efficiency based on clusters moving on circular or¬ 
bits were undoubtely over-simplified and leading to sig¬ 
nificantly overestimated values of the dynamical braking 
time-scales. This means that massive globulars in triaxial 
galaxies have probably suffered a lot of dynamical brak¬ 
ing and have reached the centre of the mother galaxy were 
they can merge to form a super-massive object (not nec¬ 
essarily a black hole) or can feed a pre-existent one. Of 
course this nucleus, if massive enough, can act in a way 
to shatter incoming globulars before they are totally or- 
bitally decayed CD examined in details the two contempo¬ 
rary effects and found that, assuming as typical globular 
cluster masses 10^, 10®, and 10^ Mq, nuclei as massive as 
5 X 10® Mq, 2 X 10® M© and 5 x 10® M© are, respectively, 
needed to effectively halt the infall of globular clusters to 
the potential minimum. 

Two equally interesting scenarios (quantitatively sup¬ 
ported in CD) are open: i) a triaxial galaxy without a pri¬ 
mordial massive nucleus can drive the merging of a signif¬ 
icant mass in the form of orbitally decayed massive globu¬ 
lars in a time scale of the order of few 10® yrs, eventually 
leading, with modes which are not trivial to be studied, 
to a central object massive enough to stop further mass 
infall; ii) a moderately massive primordial nucleus is fed 
by decayed globulars such to produce a gravitational lu¬ 
minosity in the range of normal AGNs and to grow fur- 
therly until a steady state is reached. If the mass of the 
primordial nucleus is large enough, dynamical friction on 
globular clusters is overwhelmed by the tidal shattering, 
and, moreover, the massive central object also changes the 
orbital structure around it. 

2.1. Is dynamical friction an effective cause of GCS evo¬ 
lution ? 

The actual role of dynamical friction in galaxies has been 
questioned on the basis of various arguments which actu¬ 
ally apply just to CGSs in M31 and M87, two of the best 
studied cases. It is worth nothing that even if it is plau¬ 
sible that its role is important whenever clusters move 
on almost radial orbits, the quantitative definition of the 
relevance of dynamical friction and tidal destruction has 
been done just in a triaxial galaxy (Pesce et al, 1992; CD; 
□striker et al, 1989). 

Now, while it is now quite accepted that the inner part 
of M31 is triaxial (Lindblad 1956; Stark 1977; Bertola et 


al. 1991) and so it is a galaxy where dynamical friction and 
tidal distruption effects should be significantly enhanced, a 
triaxiality for M87 is not evident. Good CCD photometric 
data by Zeilinger et al. (1993) for the inner M87 region 
show almost round isophotes in the inner M87 region (r < 
3") and a twisting occurs at 3" from the centre, the major 
axis being shifted to a position perpendicular to the jet 
and the ellipticity grows up to 0.2 at r ~ 80". This means 
that M87 is not necessarily one of the best candidate to 
investigate about the evolution of CCS distribution. 

With regard to M31 two serious observational points are: 
i) the M31 galactic nucleus seems to be significantly redder 
than globular clusters (Surdin and Charikov, 1977), ii) 
M31 globular clusters seem to show a trend of increasing 
metallicity toward the galactic centre (Huchra et al 1991), 
anyway this trend , see Fig. 2 in van den Bergh (1991) is 
quantitatively questionable. 

Let us explain why in our opinion points i) and ii) 
are much less serious indication against the importance of 
dynamical friction and GCS evolution mechanism than it 
is superficially thought. 

First of all the above mentioned data do not consti¬ 
tute a significant sample to extract general conclusions, 
referring just to one galaxy, anyway a trend of redder in¬ 
tegrated colours towards the centres is a common feature 
of many galaxies (Gallagher et al. 1980), and needs in any 
case an explanation. 

According to various authors, due to the apparent high 
metallicity of central region of M31, the decayed high mass 
clusters should have been more metal abundant (redder) 
than the only presently observed, and this needs a corre¬ 
lation between mass and metallicity for globular clusters. 
This metal abundance-mass correlation for globular clus¬ 
ters is claimed to be ad hoc because of the poor correla¬ 
tion presently observed between metal content and total 
luminosity of galactic globulars. There is an important 
caveat before concluding that also a Z mass-correlation 
is not holding. It comes by the fact that the Z-L corre¬ 
lation cannot be considered exactly representative of a Z 
mass-correlation because the mass-luminosity ratio actu¬ 
ally depends on the metal content, and it increases with 
Z. 

This means that a flat Z-L correlation transforms into 
a (more or less steep) increasing M{Z). Moreover, what 
is observed now is the present luminosity-Z correlation, 
which, in the case of efficiency of the dynamical friction 
braking is biased (respect to the initial) towards lower 
luminosities (masses) and metallicities, likely hiding an 
initial stronger correlation. 

A mass-metallicity relation for globular clusters is, 
anyway, not merely an ad hoc hypothesis, being verified 
for instance, in galaxies: the brighter galaxies contain red¬ 
der, in the average, globulars (see van der Bergh 1991 for 
M31 and our Galaxy’s clusters) and has a physical inter¬ 
pretation on the basis of a steeper, with increasing mass. 
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potential well to be overcome by the enriched material 
expelled by SNs. 

Another point that seems hardly compatible with the 
claimed efficiency of dynamical friction (acting more on 
massive globulars) is that in M87 no dependence of the 
globular cluster luminosity function on galactocentric dis¬ 
tance is found. 


2.2. The radial dependence of the GCS luminosity func¬ 
tion 


There are various reasons why the lack of evidence of a ra¬ 
dial dependence of GCS luminosity function is not a signi¬ 
ficative point against dynamical friction to be occurred on 
clusters: i) even if a spatial trend of the luminosity func¬ 
tion is present (massive globulars moved to inner regions) 
it is expected to occur in quite central regions (within the 
bulge star core radius, see Capuzzo-Dolcetta and Tesseri 
1996) not easily covered by proper observation; ii) pro¬ 
jection effects weaken any radial trend of the luminosity 
function; iii) dynamical friction reduces the average galac¬ 
tocentric distance of massive GCs more than that of light 
globulars, but this contemporarily means they are shifted 
to inner galactic zones where they likely lose their individ¬ 
uality because they become hardly observable and more 
easily destroyed by the intense tidal field. 

Let us give some quantitative support to point ii). 
Suppose to have a sample of globular clusters whose 
mean mass < m > varies with the galactocentric dis¬ 
tance in a way to have smaller masses in external re¬ 
gions {< m > varies from fO® Mq to fO® M© going 
from the centre to 5 times the core radius) and with a 
mass spectrum corresponding to a gaussian V-magnitude 
function characterized by a dispersion around the mean 
magnitude {{M/L)vq = 1.6 is assumed) which is larger 
(cTy = 3) in peripheral galactic regions than around the 
centre (uy = 1.5) (this is what qualitatively expected 
when dynamical friction and tidal disruption have been 
effective). Assuming the GCS distributed spherically ac¬ 
cording to the modified Hubble profile 


n(r) = no 



( 1 ) 


we can compare the volume LF with the projected LF, 
sampled at various distances from the centre in a galaxy 
with a distance modulus (m — M)v = 31.3 (similar to 
Virgo cluster) (see Fig. 1 a,b). 

Projection should reduce the difference among the peak 
magnitudes and the widths of the LF sampled at various 
galactocentric distances. Actually, Fig. 1 shows how a V- 
peak difference of 4 mag reduces to just 2.5 mag, while 
the width of the projected LF results almost constant (in¬ 
dependent on the galactocentric distance). 

We conclude that a radial dependence of the width of such 
LFs would be undetectable, while the detection of a varia¬ 
tion in the V-peak would require globular cluster sampling 


well within the galaxy core, because (as it is seen in Fig. 
1 b) the V-peaks of the LF sampled at r = Tc and r = 5rc 
differs for a quantity similar to the standard deviation 
of the mean (ct^ ~ 0.1 mag for a typical total sample of 
« 500 clusters). This explains why, even if evolutionary 
effects have been active on GCSs, LF radial trends have 
not been detected in the past. Higher resolution observa¬ 
tions to have larger sample abundances in inner galactic 
regions are needed. 



Fig. 1. panel a: normalized V-magnitude volume luminos¬ 
ity functions of the GCS at various galactocentric distances 
{r/vc — 0,0.25,0.5,1,5); panel b: projected luminosity func¬ 
tions at the same (projected) galactocentric distances as in 
panel a. 


3. Present and initial radial distributions 

A way to estimate the number of globular clusters lost 
during the evolution of a globular cluster system was sug¬ 
gested by McLaughlin (1995) who applied it to the M87 
galaxy. This estimate is based on the assumption that the 
stellar bulge and the globular cluster distributions were 
initially the same (due to a co-eval formation) and on 
that the shape of the stellar distribution has remained 
unchanged. The first hypothesis is supported by the ob¬ 
served similiarity of the stellar bulge and globular cluster 
projected profiles in the outer regions (outside a certain 
distance r from the centre) of various galaxies. 

The second hypothesis stands firmly on that the bulge is 
a collisionless system. The initial globular cluster distri¬ 
bution no{r) (assumed to be equal in shape to the present 
stellar distribution) is, practically, obtained by a scaling of 
the present GCS distribution to the stellar one. Of course 
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also the projected initial density profiles, cro(r), are as¬ 
sumed to be the same. Once no(r) or ao(r) have been de¬ 
termined, the number of missing clusters is given by the 
integral over the whole galaxy of the difference between 
no and n (or (Jq and cr). 

3.1. The data sets and interpolations 

We consider the well established data sets of CGS in our 
Galaxy, M31 and M87, this latter mainly for the sake of 
a comparison with McLaughlin’s (1995) results. 

To fit the stellar distribution of M31 and M87 we used 
the model of de Vaucouleurs(1958) and for our Galaxy 
the Young’s model (1976); good fits to the globular clus¬ 
ter distributions are obtained by means of the empirical 
King’s models (King 1962). 

To obtain the initial globular cluster distribution we verti¬ 
cally shift the stellar distribution to match, in the external 
region (r > f), the GGS distribution, so to have a scal¬ 
ing factor, d, depending on r. Thus, we can represent the 
initial globular cluster distribution by: 

log no = logn^ -b d{f) (2) 



- 0.5 0 0.5 1 1.5 


Log r 

Fig. 2. The globular cluster iuitial distributiou and the present 
one for our Galaxy. The dashed line refers to the point where 
the vertically shape bulge and GCS profile start to overlap 


logCTo = logo's +d{r) 


4 . 2 . Andromeda (M31) 


were na{r) and o's(r) are the (observed) stellar distribution 
volume and surface densities. Hereafter the volume and 
surface number densities will be given in kpc~^ and kpc~^ 
respectively. 

4. Number of clusters missing in our Galaxy and 
M31 

4 . 1 . Our Galaxy 

The most complete set of data of galactic globular clusters 
is still given by Webbink (1985). It refers to 154 globular 
clusters. For our purposes we need only the distance of 
each cluster from the galactic centre. A good fit to the 
present distribution is: 
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1 

0.51 
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where 6 = 49. 

As distribution of bulge stars in our Galaxy we use the 
Young’s model (1960) (see Fig. 2). Thus we obtain the 
initial globular cluster distribution: 


no(r) = 426exp 


-7.669 • 


\ ® / T’ 


-0.875 


( 5 ) 


The estimated number of globular clusters lost for our 
Galaxy is TV; = 56, i.e. about 36% of the present sample’s 
abundance. 


Various compilations of data for globular clusters in M31 
are available. We refer to the ’’Adopted Best Sample” of 
Battistini et al. (1993) compilation, for it is the most com¬ 
plete and best discussed source of data for the radial dis¬ 
tribution of globular clusters in this galaxy. 

The flattening of the M31 globular cluster distribution 
compared to that of the star spheroidal component was 
first noted by de Vaucouleurs & Buta (1978), later ques¬ 
tioned by Wirth, Smarr & Bruno (1985) who claimed a 
large part of this flattening as being due to incomplete¬ 
ness. The completeness of the Battistini et al. data in the 
inner bulge region is well addressed, and a residual flat¬ 
tening of their adopted samples respect to the spheroidal 
star component is evident. 

Matching the globular clusters data with a King’s 
model leads to the following analitycal fit to the present 
globular cluster distribution: 


a{r) = 16.7 
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0.51 
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1.52 


( 6 ) 


were a = 30. The data for the stellar distribution are 
taken from de Vaucouleurs (1958), taking the luminosity 
distribution as representative of the stellar density distri¬ 
bution (see Fig. 3). We obtain the initial globular clusters 
distribution: 

cro(r) = 785.2 • exp [-7.427 • -L 3.l] (7) 
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distribution profile is evolved, i.e. the estimate of r. A rel¬ 
ative error in r induces an error in Ni{f): 

ANi _ dNi f Af 

Ni dr Ni f ’ 
being: 

= In 10 • 27r • 10'^^'’^ • [ as(r)r dr (9) 

or Jo 

or, when the spatial density is available: 


= In 10 • 47 r • 10 '^^'’^ • d'{r) [ ng{r)r‘^ dr ( 10 ) 

or Jo 

where d'{f) is the derivative of d{f) with respect to r. 

For M31 and our Galaxy we find that the error Afjf re¬ 
flects in relative errors ANiJNi given by 0.75Af/f and 
0.63Ar/f, respectively. 


Fig. 3. Initial and present globular cluster distributions in M31 


The surface integral of cro(r) — cr(r) gives Ni ci 76 as 
number of globular clusters lost; this is 25% of the present 
number. 

4-3. Evaluation of globular cluster mass fallen to the 
galactic centres 

An approximate value of the mass fallen to the centre of 
M31 and our Galaxy can be given by mean of the knowl¬ 
edge of Ni and of the average mass of destroyed globular 
cluster < mi >. The determination of < to/ > requires a 
detailed evaluation of the tidal distruption and dynamical 
friction effects on an assumed initial mass function. By 
the way, due to that the two phenomena erode the CGS 
on apposite sides of the mass function, the mean value 
of the globular cluster mass < to > is not expected to 
change very much in time whenever the initial mass func¬ 
tion is not too asymmetric, and thus it can be chosen as a 
good reference value for < to/ >. This is confirmed by re¬ 
sults obtained with a theorical model (Capuzzo-Dolcetta 
& Tesseri 1996) under the hypotesis of a constant initial 
initial mass function (<I>(to) = constant). 

The knowledge of the mean mass of globular clusters, 
< TO >= 3.2 • 10 ®Mq for our Galaxy and < to >= 2.7 ■ 
10®M© for M31, gives as mass lost M/ = 1.8 ■ IO^Mq 
and Ml = 2.1 • IO^Mq respectively. These values should 
be compared with the nucleus masses in our galaxy (3 • 
1 O®M 0 , see Krabbe et al 1995) and in M31 (lO^M©, see 
e.g. Melia 1992). 

4-4- Sources of error 

A significant source of error evaluation of Ni is the inde¬ 
termination in the of the region where the globular cluster 


5. M87 


For the sake of comparison with previous work (McLaugh¬ 
lin 1995) we applied our method to M87. The data are 
taken from McLaughlin (1995) and from de Vaucouleurs 
and Nieto (1978, 1979) for globular clusters and halo stars, 
respectively. The fits we obtained from those distributions 
are: 


air) = 15.5 
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1 

0.5 

1 

0.51 

1 

Li + (i«)L 


4 + ( 13 ) - 
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( 11 ) 


ao{r) = 67.62 • exp 


/ V \ 0.349 


( 12 ) 


where c = 60. 

We have fixed r as the point where the two distribu¬ 
tions clearly show the same shape, as it is shown in Fig. 
4 (r = 12.5fcpc). 

The number of globular clusters lost is found to be 
Ni = 4018, which is about the same number of globular 
clusters presently observed. 

We have numerically evaluated the error induced on Ni 
by an error in r, finding that an error Afjf reflects in a 
relative errors for ANifNi ~ 2.12Ar/f; i.e. much greater 
than for our Galaxy and M31. 

This large sensitivity of Ni on the choice of f is part of 
the explanation of the great difference between our value 
of Ni and that given by McLaughlin (Ni ~ 1150). This 
difference is indeed accounted by: 

i) we have considered as best value for the radius 
r ~ 12.5 kpc instead of 8 kpc (used by McLaughlin 1995): 
it seems that the two distributions have a very similar 
shape only outside 12.5 kpc. The value of 8 kpc for f 
seems an underestimate. Infact as shown by Fig. 5, the dis¬ 
tribution obtained shifting vertically the bulge star profile 





6 


The globular cluster distributions 



2 10 1 
Log r 


Fig. 4. The initial and the present globular cluster distribution 
in M87 

to intersect the present GCS distribution at f = 8 kpc rep¬ 
resents an acceptable initial distribution for the GCS just 
if we accept as realistic that in external regions (r > 8 kpc) 
there are at present more clusters than initially. This im¬ 
plies the existence of a mechanism which populated the 
external regions, while dynamical friction and tidal dis- 
truption depopulated the inner regions. 

ii) the analytical fits are obtained with different func¬ 
tions: McLaughlin used, for both globular clusters and 
star bulge isotropic, single mass King’s (1966) models, 
while we used an empirical King’s (1963) model to fit 
the present globular cluster distribution and the de Vau- 
couleur’s (1958) model (more peaked than the King’s 1966 
one) to fit the bulge distribution. Anyway the difference in 
the analytical fits can account for just a 30% difference in 
Ni. Infact, we obtain Ni ci 1590 instead of Ni ci 1150 inte¬ 
grating our distribution up to the same radius (r = 8 kpc) 
used by McLaughlin. 

Now, if we take as mean mass of the globular clusters 
in M87 < m >= 6.6 • IO^Mq (McLaughlin 1995) our esti¬ 
mate of the mass lost is 2.65 • 10 ®Mq. If we compare this 
value with the nucleus mass of M87 (which is estimated 
to be ~ 2.4 • 10®Mq within 18pc of the nucleus, see Ford 
and al 1994) we see, also in this case, that the distruption 
of globular clusters could have strongly influenced the for¬ 
mation of the central nucleus of this galaxy. 

6. Conclusions 

It is both a reasonable and simple hypothesis that the 
globular cluster and spheroidal components of galaxies 
formed contemporarily during the first stages of proto¬ 
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Fig. 5. As in Fig. 3, but with r — 8 kpc 

galaxy collapse so to have, initially, the same spatial dis¬ 
tribution. Observed (and kinematic) spatial differences 
should be explained on the basis of evolution of the glob¬ 
ular cluster system (GCS). 

We have given reference to quantitative studies which 
point out the role of dynamical causes of this evolution in 
galaxies where clusters move on sufficiently radial orbits. 
We have also explained why most of the observational data 
available is, at present, at all unsufficient to rule out that 
such an evolution occurred. To state something meaning¬ 
ful, observations of clusters in the innermost regions (i.e. 
within the bulge core) to compare with clusters in outer 
regions of their parent galaxy, as well as kinematic data 
to determine the cluster velocity ellipsoid are needed. 

Through the comparison between the globular clus¬ 
ter and spheroidal radial distributions we determined the 
number, Ni, of clusters lost in our Galaxy, M31 and M87. 
We found that the GCSs of our Galaxy, M31 and M87 
should have been initially 1.45, 1.25, and 2 times more 
popolous than now. 

The mass of missing clusters has likely gone to the centre 
of the parent galaxy, where it can contributed to enrich 
the nucleus by an amount of the order of < m >, where 
< m > is the (present) mean value of cluster mass. This 
corresponds to 1.8 • IO^Mq, 2.1 • IO^Mq and 2.7 • 10®M© 
for the Galaxy, M31 and M87, respectively. It is relevant 
noting that these values are all very similar to available 
estimates of the nucleus masses in these galaxies. 
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